December, 2002 



Hadronic D Decays Involving Scalar Mesons 

Hai-Yang Cheng 

Institute of Physics, Academia Sinica 
Taipei, Taiwan 115, Republic of China 



Abstract 

The nonleptonic weak decays of charmed mesons into a scalar meson and a pseudoscalar meson 
are studied. The scalar mesons under consideration are a [or /o(600)], n, /o(980), ao(980) and 
Kq(1430). A consistent picture provided by the data suggests that the light scalars below or near 
1 GeV form an SU(3) flavor nonet and are predominately the q 2 q 2 states, while the scalar mesons 
above 1 GeV can be described as a qq nonet. Hence, we designate q 2 q 2 to a, k, a (980), / (980) and 
qq to Kq . Sizable weak annihilation contributions induced from final-state interactions are essential 
for understanding the data. Except for the Cabibbo doubly suppressed channel D + — > /qK + , the 
data of D —* air, fair, foK, KqTt can be accommodated in the generalized factorization approach. 
However, the predicted rates for D —* aovr, clqK are too small by one to two orders of magnitude 
when compared with the preliminary measurements. Whether or not one can differentiate between 
the two-quark and four-quark pictures for the /o(980) produced in the hadronic charm decays 
depends on the isoscalar fo — a mixing angle in the qq model. 
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I. INTRODUCTION 



There are two essential ingredients for understanding the hadronic decays of charmed 
mesons. First, the nonfactorizable contributions to the internal ^-emission amplitude, which 
is naively expected to be color suppressed, is very sizable. Second, final-state interactions 
(FSIs) play an essential role. The nonfactorizable corrections to nonleptonic charm decays 
will compensate the Fierz terms appearing in the factorization approach to render the naively 
color-suppressed modes no longer color suppressed. The weak annihilation (W -exchange 
or ^-annihilation) amplitude receives long-distance contributions via inelastic final-state 
interactions from the leading tree or color-suppressed amplitude. As a consequence, weak 
annihilation has a sizable magnitude comparable to the color-suppressed internal lU-emission 
with a large phase relative to the tree amplitude. A well known example is the decay 
D° — > K (f) which proceeds only through the VF-exchange process. Even in the absence of 
the short-distance W^-exchange contribution, rescattering effects required by unitarity can 
produce this reaction [1]. Then it was shown in [2] that this rescattering diagram belongs to 
the generic lU-exchange topology. 

There exist several different forms of FSIs: elastic scattering and inelastic scattering 
such as quark exchange, resonance formation,- • •, etc. The resonance formation of FSIs 
via qq resonances is probably the most important one to hadronic charm decays owing to 
the existence of an abundant spectrum of resonances known to exist at energies close to 
the mass of the charmed meson. Since FSIs are nonperturbative in nature, in principle 
it is notoriously difficult to calculate their effects. Nevertheless, most of the properties of 
resonances follow from unitarity alone, without regard to the dynamical mechanism that 
produces the resonance [3,4]. Consequently, the effect of resonance-induced FSIs can be 
described in a model-independent manner in terms of the masses and decay widths of the 
nearby resonances (for details, see e.g. [5]). 

It has been established sometime ago that a least model-independent analysis of heavy 
meson decays can be carried out in the so-called quark-diagram approach [6,7,2]. Based on 
SU(3) flavor symmetry, this model-independent analysis enables us to extract the topological 
quark-graph amplitudes for D — > PP, VP decays and see the relative importance of different 
underlying decay mechanisms. From this analysis one can learn the importance of the weak 
annihilation amplitude and a nontrivial phase between tree and color-suppressed amplitudes 
[8]. 

In the present paper we shall study the nonleptonic decays of charmed mesons into a 
pseudoscalar meson and a scalar meson. Light scalar mesons are traditionally studied in low 
energy S'-wave 7T7T, Kn and KK scattering experiments and in pp and ZN annihilations. 
Thanks to the powerful Dalitz plot analysis technique, many scalar meson production mea- 
surements in charm decays are now available from the dedicated experiments conducted at 
CLEO, E791, FOCUS, and BaBar. Hence the study of three-body decays of charmed mesons 
opens a new avenue to the understanding of the light scalar meson spectroscopy. Specifi- 
cally, the decays D — > f 7i(K), D — > a ir(K), D — > K* rr and D + — > <jit + have been observed. 
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Moreover, in some three-body decays of charmed mesons, the intermediate scalar meson ac- 
counts for the main contribution to the total decay rate. For example, — > /o(980)7r + and 
P+ — > / (1370)7T + account for almost 90% of the P+ — > rate [9], while about half 

of the total decay rate of D + — > 7r + 7r + 7r~ comes from D + — > cnr + [9]. 

The study of D — > SP is very similar to P — > PP except for the fact that the quark 
structure of the scalar mesons, especially / (980) and a (980), is still not clear. A consistent 
picture provided by the data implies that light scalar mesons below or near 1 GeV can be 
described by the q 2 q 2 states, while scalars above 1 GeV will form a conventional qq nonet. 
Another salient feature is that the decay constant of the scalar meson is either zero or very 
small. We shall see later that final-state interactions are essential for understanding the 
D — > SP data. It is hoped that through the study of D — > SP, old puzzles related to the 
internal structure and related parameters, e.g. the masses and widths, of light scalar mesons 
can receive new understanding. 

This work is organized as follows. In Sec. II we summarize the experimental measure- 
ments of D — > SP decays and emphasize that many results are still preliminary. We then 
discuss the various properties of the scalar mesons in Sec. Ill, for example, the quark struc- 
ture, the decay constants and the form factors. Sec. IV is devoted to the quark-diagram 
scheme and its implication for final-state interactions. We analyze the D — > SP data in Sec. 
V based on the generalized factorization approach in conjunction with FSIs. Conclusion is 
made in Sec. VI. 



II. EXPERIMENTAL STATUS 

It is known that three-body decays of heavy mesons provide a rich laboratory for study- 
ing the intermediate state resonances. The Dalitz plot analysis is a very useful technique 
for this purpose. We are interested in D — > SP (S: scalar meson, P: pseudoscalar me- 
son) decays extracted from the three-body decays of charmed mesons. Some recent re- 
sults (many being preliminary) are available from E791 [9], CLEO [10], FOCUS [11] and 
BaBar [12]. The + scalar mesons that have been studied in charm decays are <r(500) [or 
/o(600)], / (980), /o(1370), a (980), a (1450), re and K*(U30). The results of various exper- 
iments are summarized in Table I where the product of B(D — > SP) and B(S — > P1P2) 
is listed. In order to extract the branching ratios for D — > foP, we use the value of 
r(/ -> 7T7r)/[r(/ -> tttt) + r(/ -> KK)\ = 0.68 [13]. Therefore, 

B(/o(980) -> K+K-) = 0.16, B(f (980) -> tt+tt-) = 0.45 . (2.1) 

For D — > a P, we apply the PDG (Particle Data Group) average, T(a — > KK)/T(a — > 
nr]) = 0.177 ± 0.024 [14], to obtain 

£(a+(980) -> K + K°) = 8(^(980) -> K~K°) = 0.15 ± 0.02 , 

B(a° (980) -> K + K-) = 0.075 ± 0.010 . (2.2) 
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Needless to say, it is of great importance to have more precise measurements of the branching 
fractions of fo and do- 

TABLE I. Experimental branching ratios of various D — > SP decays measured by ARGUS, 
E687, E691, E791, CLEO, FOCUS and BaBar, where use of Eqs. (2.1) and (2.2) for the branching 
fractions of /o(980) and ao(980) has been made. For simplicity and convenience, we have dropped 
the mass identification for / (980), a (980) and if *(1430). 
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Several remarks are in order. 

1. The Cabibbo doubly suppressed mode D + — > K + K + K~ has been recently observed 
by FOCUS [11]. The Cabibbo-allowed D° — > 0^7^^ and doubly Cabibbo-suppressed 
mode -D — > agi^ + have been extracted from the three-body decay D° — > K + K~K° 
by BaBar [12]. 

2. The decay D + — > A'~7r + 7r + has been measured by E691 [15] and E687 [16] and the 
combined branching ratio for D + — > Kq°tt + is quoted to be (3.7 ± 0.4)% by PDG [14] 
(see also Table I). A highly unusual feature is that this three-body decay is dominated 
by the nonresonant contribution at 90% level, whereas it is known that nonresonant 
effects account for at most 10% in other three-body decay modes of charmed mesons 
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[14,5]. A recent Dalitz plot analysis by E791 [9] reveals that a best fit to the data is 
obtained if the presence of an additional scalar resonance called k is included. As a 
consequence, the nonresonant decay fraction drops from 90% to (13 ± 6)%, whereas 
ktt + accounts for (48 ± 12)% of the total rate. Therefore, the branching ratio of 
D+ -> K* °n+ is dropped from (3.7 ± 0.4)% to (1.8 ± 0.3)%. We shall see in Sec. V 
that the form factor for D — > Kq transition extracted from the E791 experiment is 
more close to the theoretical expectation than that inferred from E691 and E687. 

3. The new CLEO and BaBar results on the Cab ibbo- allowed decay D° — > fo~K° are 
consistent with the early measurements by ARGUS [17] and E687 [16] quoted in Table 
I from PDG. The Cabibbo doubly suppressed mode D + — > foK + was first measured 
by FOCUS recently. 

4. There are four measurements of D° — > Kq~tt + : three from D° — > K°7r + 7T~ by ARGUS 
[17], E687 [16], CLEO [10], and one from D° -> K~ir + ir° by CLEO. The CLEO result 
(1.17±0.26)% for the branching ratio of D° — > Kq~it + extracted from D° — > K~7t + tt° 
is in good agreement with ARGUS and E687 (see Table I), while the CLEO number 
(7.0l?;a) x 10^ 3 determined from D° — > X°7r + 7r _ is slightly lower. 

5. As for — > /o 7r+ ) four measured results by E687, E791 and FOCUS are shown in 
Table I. The old measurement by E687 and two new ones by FOCUS are larger than the 
E791 one. The preliminary FOCUS measurement indicates that the /o(980) resonance 
accounts for (94.4 ±3.8)% of the total — > it + it + it~ rate [11]. Later we shall use the 
average value B(D+ -> f n + ) = (1.8 ± 0.3)% in Table IV. 

6. As stressed in the Introduction, there exist three-body decay modes that are domi- 
nated by the scalar resonances. Apart from the decays Df — > fo7r + and D + — > <jtt + 
as mentioned in the Introduction, some other examples are — > fo(980)K + and 
D + — > / (980)-ft' + which account for 72% and 44.5%, respectively, of the decays 
Df -> K + K+K~ and D + -> K+K+R- [11]. 

7. The production of the resonance /o(1370) in D° — > K°7t + tt~ — > / (1370)!K , D + — > 
K + K~7i + — > /o(1370)7r + and — > 7r + 7r + 7r _ — > /o(1370)7r + has been measured by 
ARGUS, E687, CLEO, by FOCUS and by E791, respectively. Since the branching 
fractions of / (1370) into n + 7i~ , K + K~ are unknown, we will not discuss it until Sec. 
V.D. 

8. Some preliminary measurements of D — > SP do not have yet enough statistical signif- 
icance, for example, the decays D — > X + and D+ -> K K + . 
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III. PHYSICAL PROPERTIES OF SCALAR MESONS 



The masses and widths of the + scalar mesons relevant for our purposes are summarized 
in Table II. The a meson observed in D + — > 7r + 7r + 7r~ decay by E791 [9] has a mass of 
478l|| ± 17 MeV and a width of 324±|g ± 21 MeV. Recently, the decay D° -> fT°7r+7r- 
has been analyzed by CLEO [10]. By replacing the nonresonant contribution with a Kga 
component, it is found that m a = 513 ± 32 MeV and T a = 335 ± 67 MeV [10], in accordance 
with E791. The isodoublet scalar resonance k observed in the decay D + — > K~it + -k + by 
E791 has a mass of 797 ± 19 ± 43 MeV and a width of 410 ± 43 ± 87 MeV [9]. However, the 
signal of k is much less evident than a. Indeed, this resonance is not confirmed by CLEO 
in the Dalitz analysis of the decay D° — > K^-k + -k° [10]. The well established scalars /o(980) 
and a (980) are narrow, while a and k are very broad. 

TABLE II. The masses and widths of the 1 3 Pq scalar mesons (except for k, see the mini-review 
in [18]) quoted in [14]. 





a 




K 


/o(980) 


a (980) 


^(1430) 


mass 400 — 


1200 MeV 


700- 


900 MeV 


980 ± 10 MeV 


984.7 ± 1.2 MeV 


1412 ± 6 MeV 


width 600 - 


1000 MeV 


400- 


600 MeV 


40- 100 MeV 


50- 100 MeV 


294 ± 23 MeV 



A. Quark structure of scalar mesons 

It is known that the identification of scalar mesons is difficult experimentally and the 
underlying structure of scalar mesons is not well established theoretically (for a review, see 
e.g. [18-20]). It has been suggested that the light scalars-the isoscalars <r(500), / (980), 
the isodoublet k and the isovector a (980)-form an SU(3) flavor nonet. In the naive quark 
model, the flavor wave functions of these scalars read 

fo = ss, 

ciq = ud, <2q = du, (3-1) 
R° = sd, k~ = su. 

However, this model immediately faces two difficulties:* (i) It is impossible to understand 
the mass degeneracy of / (980) and a (980). (ii) It is hard to explain why a and k are 



*However, for a different point of view of these difficulties in the qq picture, see e.g. [21]. 



a = —j={uu + dd), 
a/2 



a 



o 



1 



(uu — dd), 



°~ a/2 v 

K + = US, K° = ds, 
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broader than / (980) and a (980). Recalling that a — > ni], a — > nn and k — > fT7r are OZI 
allowed (but not OZI superallowed !) while /o — > 7r7r is OZI suppressed as it is mediated by 
the exchange of two gluons, it is thus expected that m K ^> T K ~ ~ T ao > r^ , a relation 
not borne out by experiment. 

Although the data of D+ -> /o(980)tt + and -> / (980)7 imply the copious / (980) 
production via its ss component, there are some experimental evidences indicating that 
/o(980) is not purely an ss state. First, the measurements of J/ip — * / O (98O)0 and J/ip — > 
/o(980)u; 

B(J/ip -> / O (98O)0) = (3.2 ± 0.9) x 10" 4 , 

B(J/ip -> / (980)w) = (1.4 ± 0.5) x 10" 4 (3.2) 

clearly indicate the existence of the non-strange and strange quark content in / (980). Sec- 
ond, the fact that / (980) and a (980) have similar widths and that the f width is dominated 
by nn also suggests the composition of uu and dd pairs in / (980); that is, / (980) — > nn 
should not be OZI suppressed relative to a (980) — > nrj. Therefore, isoscalars a and f must 
have a mixing 

fo = ss cos 9 + nn sin 9, a = —ss sin 9 + nn cos 9, (3.3) 

with nn = (uu + dd) / y/2. 

The a — /o(980) mixing angle can be inferred from the decays J/ip —> /o(98O)0 and 
J/tP -> / (980)u; [22] t 

B( W - /o(980V) 1 2 



B{JH> - /o(98O)0) A 



tan 2 #, (3.4) 



where the deviation of the parameter A from unity characterizes the suppression of the ss 
pair production; that is, A = 1 in the SU(3) limit. From the data (3.2) we obtain 

#=(34 ±6)°, or # = (146 ±6)° (3.5) 

for A = 1. Another information on the mixing angle can be obtained from the /o(980) 
coupling to nn and KK [22] : 



R g = 9 \ K+K ~ = I(A + v^cot#) 2 . (3.6) 



Tt has been shown by the / (980) production data in Z° decays at OPAL [23] and DELPHI [24] 
that /o(980) is composed essentially of uu and dd pairs. This favors a mixing angle close to n/2. 
However, it is in contradiction to the experimental observation that the final state /o (980)0 in 
hadronic J/ip decays has a larger rate than /o(980)oj. 
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Using the average value of R g = 4.03 ± 0.14 obtained from the measurements: 4.00 ± 0.14 
by KLOE [25], 4.6 ± 0.8 by SND [26] and 6.1 ± 2.0 from CMD-2 [27], we find 

9 = (25.1 ±0.5)°, or 9 = (164.3 ±0.2)° (3.7) 

for A = 1. However, the WA102 experiment on / (980) production in central pp collisions 
yields a result R g = 1.63 ± 0.46 [28], which differs from the aforementioned three measure- 
ments. This leads to 

= (42.3ljj;ij) , or = (158±2)°, (3.8) 

again for A = 1. 

Recently, a phenomenological analysis of the radiative decays — > f (980)7 an d 
/o(980) -> 77 yields 

= (5±5)°, or 0=(138±6)° (3.9) 

with the second solution being more preferable [29]. In this analysis, — > fo (980)7 is 
calculated at the quark level by considering the ss quark loop coupled to both and /o(980).* 
However, the experimental analysis and the theoretical study of this radiative decay are 
practically based on the chiral-loop picture, namely, — > K + K~ — > K + K~ r y — > f (980)7. 
It turns out that the predicted branching ratio in the qq picture is at most of order 5 x 10~ 5 
(see e.g. [31]), while experimentally [14] 

5(0 - / (980) 7 ) = (3.318;|) x 10~ 4 . (3.10) 

This is because the /o(980) coupling to K + K~ is not strong enough as in the case of the 
four-quark model to be discussed shortly where fo — > K + K~ is OZI superallowed. 

In short, it is not clear if there exists a universal mixing angle which fits simultaneously 
to all the measurements from hadronic J/vp decays, the /o(980) coupling to 7t + tt~ and K + K~ 
and the radiative decay — > fo (980)7 followed by fo — > 77. § Eqs. (3.5) and (3.9) indicate 



*It is pointed out in [30] that the mechanism m ss —> ssj — ► /o (980)7 without creation of 
and annihilation of an additional uu pair cannot explain the /q(980) spectrum observed in — > 
7/o(980) — > 77r°7r° process because it does not contain the K + K~ intermediate state. For a 
criticism of [29], see also the remark in the footnote [28] in the first paper of [30]. 

§The analysis of the three-body decays of D+ -► /o(980)tt + , i?o°(1430)7r+ and D+ / (980)vr+ 
gives 6 = (42.I4I73) in [32]. However, in this analysis, the weak annihilation contribution has 
been neglected and SU(3) symmetry has been applied to relate foir to KqTt, a procedure which is 
not justified since fo and Kq belong to different SU(3) flavor nonets (see the discussion in Sec. IV). 
If the mass parameters m n n and m S s are assigned to the nn and ss states, respectively, one will 
have the mass relations: m^ n = m} a cos 6 2 + nij o sin 9 2 and = m 2 sin 9 2 + m 2 o cos 9 2 from Eq. 
(3.3). Inserting the dynamically generated NJL-type masses for m nn and m sS , it is found in [33] 
that 9 = ±(18 ± 2)° provided that m c = 600 MeV. 
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that ~ 140° is preferred, while 9 ~ 34° is also allowed provided that R g is of order 2. At 
any rate, the above two possible allowed angles imply the dominance of ss in the /o(980) 
wave function and nn in a. However, as stressed before, the 2-quark picture for / (980) has 
the difficulty of explaining the absolute — > fo (980)7 rate. 

As for ao(980), it appears at first sight that one needs an ss content in ao in order to 
explain the radiative decay — > ao (980)7; otherwise, it is OZI suppressed. However, since 
ao is an isovector while ss is isoscalar, the mixing of (uu — dd) with ss is not allowed in the 
a wave function within the 2-quark description.** Nevertheless, — > a (980)7 can proceed 
through the process K + K~ — > K + K~^y — * a (980)7 as both and a (980) couple to 
K + K~ . Indeed, it has been suggested that both ao(980) and /o(980) can be interpreted as 
a KK molecular bound state which is treated as an extended object. Since both /o(980) 
and a (980) couple strongly to KK as they are just below the KK threshold, they can 
be imagined as an ss and nn core states, respectively, surrounded by a virtual KK cloud 
[34]. In this KK molecular picture, one can explain the decay fo — > tttt without the light 
non-strange quark content in / (980) and the decay — > ao (980)7 without the need of an 
intrinsic strange quark component in a ; both decays are allowed by the OZI rule in the 
sense that only one gluon exchange is needed. 

However, there are several difficulties with this KK molecular picture. First, the KK 
molecular width is less than its binding energy of order 20 MeV [34], while the measured 
widths of / (980) and a (980) lie in the range of 40 to 100 MeV [14]. Second, it is expected 
in this model that £>(0 — > /o (980)7) /£>(0 ~^ a o (980)7) ~ 1> while this ratio is measured to 
be 3.8 ±1.0 [14]. (The most recent result is 6.1 ± 0.6 by KLOE [35].) Third, the predicted 
branching ratios for both — * / (980)7 and — > a (980)7 are only of order 10~ 5 [31] which 
are too small compared to (3.10) and [36] 

B((j) -> a (980)7) = (0.88 ± 0.17) x 10" 4 . (3.11) 

Alternatively, the aforementioned difficulties''"'" with ao and fo can be circumvented in the 
four-quark model in which one writes symbolically [39] 

a = udud, fo = (usus + dsds) / V2, 



**Even in the presence of the hidden ss content in ao(980) within the 4-quark model, the direct 
radiative decay — ► ao (980)7 i s prohibited owing to the opposite sign between the uu and dd 
components in ao(980). This means that it is necessary to consider the contribution from the 
K + K~ intermediate states. 

^Likewise, it has been argued in the literature that a is not a qq state [37]. Furthermore, the 
QCD sum rule calculation also indicates that the lightest scalars are nearly decoupled from qq, 
suggesting a non-qq structure [38]. In short, one always has some troubles when the light scalar 
mesons are identified as qq states. 
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a[J = (usus — dsds)/V%, 
k + = udds, k° = udus, 



Oq = usds, 
R° = usud, 



a = dsns, 
k~ = dsud. 



(3.12) 



This is supported by a recent lattice calculation [40]. This q 2 q 2 scenario for light scalars has 
several major advantages: (i) The mass degeneracy of / (980) and a (980) is natural and the 
mass hierarchy pattern of the SU(3) nonet is understandable, (ii) Why a and k are broader 
than /o and a can be explained. The decays a — > rnv, k — > Kn and / , ao — *■ KK are OZI 
superallowed without the need of any gluon exchange, while f — > nn and a — > nr] are OZI 
allowed as it is mediated by one gluon exchange. Since / (980) and a (980) are very close to 
the KK threshold, the / (980) width is dominated by the nn state and a governed by the 
ni] state. Consequently, their widths are narrower than a and k. (iii) It predicts the relation 



which is in good agreement with the experimental value of 0.44 ±0.20 [14]. (iv) The coupling 
of f and a to KK is strong enough as the strong decays f — > KK and a — > KK are 
OZI supperallowed. Consequently, the branching ratio of — > (/ , a )7 can be as large 
as of order 10~ 4 [3l]M (v) It is concluded in [30] that production of / (980) and a (980) 
in the — > 7/0 (980) — > 77r°7r° and — > 700 (980) — > 77r°?7 decays is caused by the four- 
quark transitions, resulting in strong restrictions on the large- N c expansions of the decay 
amplitudes. The analysis shows that these constraints give new evidences in favor of the 
four-quark picture of /o(980) and ao(980) mesons. 

Therefore, it appears that the four-quark state in core with the KK in the outer regime 
gives a more realistic description of the light scalar mesons. If scalar mesons near and below 
1 GeV are non-gg states, then the + mesons in the 1.3 — 1.7 GeV mass region may be 
more conventional. For example, it is natural to assume that / (1370), a (1450), Kq(1A30) 
and /o(1500)//o(1710) are in the same SU(3) flavor nonet in the states nn, ud, us and ss, 
respectively [18]. In other words, they may have a simple qq interpretation. A global picture 
emerged from above discussions is as follows: The scalar meson states above 1 GeV form a 
qq nonet with some possible mixing with glueballs, whereas the light scalar mesons below or 
near 1 GeV form predominately a qqqq nonet with a possible mixing with + qq and glueball 
states (see also [20]). This is understandable because in the qq quark model, the + meson 



"Just as in the KK molecular model, the ratio r = B(4> — > /07) — > 007) is also an issue 
in the four-quark model in which gj K+K- = 9a K+K~ an d hence B(<f) — ► /07) = B(4> — > 007) is 
predicted, in disagreement with the observed value of 3.8 ± 1.0 [14]. Close and Kirk [41] proposed 
that r can be explained by considering a large ao — /o mixing. However, as pointed out in [42], 
the isospin- violating ao — /o mixing is small, analogous to the smallness of 7r° — r\ — r( mixing, and 
its correction to r amounts to at most a few percent. One possibility for a large r is that the 
superposition of the 4-quark state and KK has a different weight in /o(980) and ao(980). 



gW -+ /o(980V) _ 1 
B(J/*P - / O (98O)0) 2' 



(3.13) 
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has a unit of orbital angular momentum and hence it should have a higher mass above 1 
GeV. On the contrary, four quarks q 2 q 2 can form a + meson without introducing a unit of 
orbital angular momentum. Moreover, color and spin dependent interactions favor a flavor 
nonet configuration with attraction between the qq and qq pairs. Therefore, the + qqqq 
nonet has a mass near or below 1 GeV. 

It is conceivable that the two-quark and four-quark descriptions of light scalars, especially 
/o(980) and a (980), may lead to some different implications for the hadronic weak decays 
of charmed mesons into the final state containing a scalar meson. This will be explored in 
Sec. V. 

B. Decay constants 

The scalar mesons under consideration are <r(500), k, / o (980), a (980) and _K"q(1430). 
The decay constants of scalar and pseudoscalar mesons are defined by 

(0\A,\P(q)) = if P q„ (0\V,\S(q)) = f s q». (3.14) 

For the neutral scalars a, f and a[j, the decay constant must be zero owing to charge 
conjugation invariance or conservation of vector current: 

U = fh = M=0- (3-15) 

Applying the equation of motion, it is easily seen that the decay constant of Kq + (ciq ) is 
proportional to the mass difference between the constituent s (d) and u quarks. Contrary 
to the case of pseudoscalar mesons, the decay constant of the scalar meson vanishes in the 
SU(3) limit or even in the isospin limit. Therefore, the decay constant of 1^(1430) and the 
charged a (980) is suppressed. We shall use the values 

f a ± = 1.1 MeV, f K . = 42 MeV (3.16) 

obtained from the finite-energy sum rules [43] . (A different calculation of the scalar meson 
decay constants based on the generalized NLJ model is given in [44].) Since they are derived 
using the qq quark model, it is not clear if the czq decay constant remains the same in the 
q 2 q 2 picture, though it is generally expected that the decay constant is suppressed in the 
latter scenario because a four-quark state is larger than a two-quark state [43]. 

As for the decay constant of n, we apply the equation of motion to Eq. (3.14) to obtain 

m lja = i(m d - m u )(a^\du\0), m 2 J K = i(m s - m u ) (k~ \su\0) , (3.17) 

and assume (k~\su\0) ks (oq | Ju|0). It follows that f K xs 65 MeV for m u = 4.8 MeV, m d = 8.7 
MeV, m s = 164 MeV [45] and m K = 800 MeV. In short, the decay constants of scalar mesons 
are either zero or very small. 



11 



C. Form factors 



Form factors for D ^ P and D — > S transitions are defined by [46] 

(p(p)\v»\D(p D )) = (p D , + P , -H^± q ^j F P p ( q 2 ) + ^^^ D V), 

{ PD, +P»- J FP S {?) + ^-^q, F DS (q 2 ) 



(S(p)\A,\D(p D )) 



, (3.18) 



where = (po — p)n- The form factors relevant for D — > SP decays are FP p (q 2 ) and 
FP s (q 2 ). For P — > P form factors, we will use the Melikhov-Stech (MS) model [47] based 
on the constituent quark picture. Other form factor models give similar results. 

As discussed in Sec. III. A, the light scalar mesons a, k, /q(980) and ao(980) are predomi- 
nately q 2 q 2 , while K£ is described by the qq state. Nevertheless, it is useful to see what are the 
predictions of D — > S form factors in the conventional quark model. There are some existing 
calculations of the form factor F® s (0) in the literature (see Table III). Paver and Riazuddin 
[48] obtained P Da (0) = 0.74(/ D /200 MeV). Gatto et al. got FP a (0) = 0.57 ±0.09 [49] using 
the constituent quark model (CQM). Based on the same model, Deandrea et al. [50] obtained 
F Ds+/o (0) = 0.641&03 assuming a pure ss state for f . A value of FP tfo (ml) = 0.36±g;8| is 
obtained by Gourdin, Keum and Pham [51] based on a fit to the old data of P+ — > / (980)7r + . 
The B — > do form factor is estimated by Chernyak [52] to be P(f ao (0) ~ 0.46. Using the 
scaling law, it leads to P Da °(0) « P Ba °(0) ^m B /m D = 0.77. Since the conventional quark 
model is not applicable to light scalars with four-quark content, we shall use the measured 
decay rates to extract the D — > S form factors (except for D — > a ) in Sec. V and the results 
are summarized in Table III. 

TABLE III. The D -> S transition form factors FP s (0) at q 2 = in various models. Except 
for the D — > ao form factor, the other form factors in this work are obtained by a fit to the data. 
The D — > /o form factor is obtained from the P+ — > /o one via Eq. (3.20). 

Transition [48] [49] [50] This work 

D^a 0.74 0.57 ±0.09 " ' 0.42 ± 0.05 

D -»• / 0.26 ± 0.02 

£>+ - /o 0.64ig;g| 0.52 ± 0.04 

D -► 0.77 
D ^ k 0.85 ±0.10 

D,D+^K$ 1.20 ±0.07 



In the qq description of /o(980), it follows from Eq. (3.3) that 

Fp« = ^sin9F D " f °\ F D+h = ±= sin 9 F^ & , F D ^ = cos£P D ^, (3.19) 
v 2 v 2 
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(a) (b) 

FIG. 1. The Df -> / (980) and D + -» /o(980) transition form factors, where / (980) is de- 
scribed by a q 2 q 2 state. 

where the superscript qq denotes the quark content of fo involved in the transition. In the 

7~)0 fuu /~)"r" fdd J~~)~^~ f s ^ 

limit of SU(3) symmetry, F /o = F /o = F s /o and hence 

Fo D 0/0 = ^ +/ ° = 4= ^o Ds+/ ° tan 5. (3.20) 
V2 

In the four-quark picture, one has (see Fig. 1) 

^ /o (0) = ^o Ds+/ °(0), (3.21) 

where use of the /o(980) flavor function ss(uu + dd) / v2 has been made. For 9 ~ 34°, we see 
that the relation between F Df ° and F 3 To is very similar in the qq and q 2 q 2 pictures, but if 
9 ~ 140° then F will have an opposite sign to F 3 in the former model. 

To compute the D —>■ S form factors using the qq model, it is worth mentioning the 
Isgur-Scora-Grinstein-Wise (ISGW) model [53] and its improved version the ISGW2 model 
[54]. Contrary to the above-mentioned form-factor calculations shown in Table III, the form 
factors of interest in the ISGW model are F® s evaluated at q 2 = q 2 ^ = {mo — ms) 2 , the 
maximum momentum transfer, recalling that Fi(0) = -Fo(O). The reason for considering 
the form factor at zero recoil is that the form-factor q 2 dependence in the ISGW model is 
proportional to Exp[— (q^ — q 2 )}. Hence, the form factor decreases exponentially as a function 
of (g^ — q 2 ). Consequently, the form factor is unreasonably small at q 2 = 0. This has been 
improved in the ISGW2 model in which the form factor has a more realistic behavior at 
large (q 2 m — q 2 ). However, we find that form factors in the ISGW2 model calculated even 
at zero recoil are already small compared to the other model calculations at q 2 = 0. For 
example, F? a ° and F x 8 at zero recoil are found to be 0.52 and 0.29 respectively in the 
ISGW model and 0.12 as well as 0.09 respectively in the ISGW2 model. These form factors 
become even much smaller at q 2 = 0. Therefore, the ISGW model and especially the ISGW2 
version predict much smaller D — > S form factors than other quark models. 

For the q 2 dependence of the D — ► S form factors, we shall assume the pole dominance: 

F DS (q 2 ) = - F \f 2 , (3.22) 
1 — q 2 /m* 
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with m* being the mass of the 0~ pole state with the same quark content of the current 
under consideration. 

In the MS [47] or the BSW [46] model, the typical D — ^ P form factors have the values 
F D7T± (0) = 0.69 and F® K (0) = 0.78. In general, it is conceivable that the form factor for 
D — > a transition is comparable to the D — > 7T° one. The argument goes as follows. If the 
scalar meson is made from qq, its distribution amplitude has the form [45] 



<Ps{x) = 6f s x(l - x) 



(3.23) 



where fs is the decay constant of the scalar meson S, B n are constants, and CfJ 2 is the 
Gegenbauer polynomial. For the isosinglet scalar mesons a and /o, their decay constants 
vanish, but the combination fsB n can be nonzero. For a, f and Oq, charge conjugation 
invariance implies that (f>(x) = —0(1 — x); that is, the distribution amplitude vanishes at 
x — 1/2. Using |/s.Bi| ~ 75 MeV obtained in [45], we have 

(f) a {x) w6Sia:(l-a:)(3-6a:), (3.24) 

where B 1 = —fsB 1 . It is clear that the a distribution amplitude peaks at x — 0.25 and 
0.75. Now the D meson wave function is peaked at x ~ A/ mo ~ 1/3 [55]. Recall that the 
asymptotic pion distribution amplitude has the familiar expression 

fa(x) = 6f v x(l-x), (3.25) 

which has a peak at x — 1/2. Though B\ is smaller than /„., it is anticipated that the D — > a 
transition form factor is similar to that of D — > 7r° one because the peak of <p a is close to 
that of the D distribution amplitude. However, it is not clear if this argument still holds for 
the scalar mesons which are bound states of q 2 q 2 . 
For _K"q(1430), the distribution amplitude reads 

(/> K5 (x) » 6f K *x(l -x)[l + fii(3 - 6x)} . (3.26) 

It is easy to check that 4>k* has a large peak at x — 0.25 and a small peak at x — 0.75. 
Consequently, it is natural to have a large D — > Kq transition form factor as shown in Table 
III. Another argument favoring a large D — > Kq form factor is as follows (see [56]). Consider 
the decay D° — > Kq~7t + and apply PCAC to evaluate the matrix element q^(KQ |s7 M (l — 
75)c|D°). Assuming that this two-body matrix element is saturated by the Df pole, we find 

2 

(n + \(ud)\0){K* o -\(sc)\D°) = iV2Uf Ds g K; - D0Dt J?* ^ (3.27) 



Next apply the SU(4) symmetry to relate strong coupling of Kq D°D+ to Kq n + K 

9k* -d°d+ = (3-28) 
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The coupling g K „^^ + ^o can be determined from the measured decay rate of Kq — > K n 
via 

where p c is the cm. momentum in the Kq rest frame. Using the Kq width given in Table 
II, we obtain g„»_ +T7 o = 4.9 GeV. Since 

(n + \(ud)\0)(K*~\(sc)\D°) = Mm 2 D -m 2 Kt) )F DK5 (ml), (3.30) 

it follows that 



r- 3k*~tt+k°^d+ m 2 n ( f n + \ 



m D - m D s - m £ \ 270 MeV , 

This is consistent with the value of 1.20 ± 0.07 extracted directly from D + — > K n + (see 



Sec. V.A). 



IV. QUARK DIAGRAM SCHEME AND FINAL-STATE INTERACTIONS 

It has been established sometime ago that a least model-independent analysis of heavy 
meson decays can be carried out in the so-called quark-diagram approach. In this diagram- 
matic scenario, all two-body nonleptonic weak decays of heavy mesons can be expressed 
in terms of six distinct quark diagrams [6,7,2]: T, the color-allowed external lU-emission 
tree diagram; C, the color-suppressed internal IF-emission diagram; E, the IF-exchange 
diagram; A, the lU-annihilation diagram; P, the horizontal WMoop diagram; and V, the 
vertical WMoop diagram. (The one-gluon exchange approximation of the P graph is the 
so-called "penguin diagram".) It should be stressed that these quark diagrams are classified 
according to the topologies of weak interactions with all strong interaction effects included 
and hence they are not Feynman graphs. All quark graphs used in this approach are topo- 
logical and meant to have all the strong interactions included, i.e. gluon lines are included 
in all possible ways. Therefore, topological graphs can provide information on final-state 
interactions (FSIs). 

Based on the SU(3) flavor symmetry, this model-independent analysis enables us to 
extract the topological quark-graph amplitudes and see the relative importance of different 
underlying decay mechanisms. The quark-diagram scheme, in addition to be helpful in 
organizing the theoretical calculations, can be used to analyze the experimental data directly. 
When enough measurements are made with sufficient accuracy, we can find out the values of 
each quark-diagram amplitude from experiment and compare to theoretical results, especially 
checking whether there are any final-state interactions or whether the weak annihilations can 
be ignored as often asserted in the literature. 
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For charmed meson decays, the penguin contributions are negligible owing to the good 
approximation V u< jV* d pz —V us V* s and the smallness of V u bV* b . Hence, for D — > PP, VP, VV 
decays, only T, C, E and A contribute. The reduced quark-graph amplitudes T, C, E, A have 
been extracted from Cabibbo-allowed D — > PP decays by Rosner [8,57] with the results: 

T = (2.67 ±0.20) x 10~ 6 GeV, 

C = (2.03 ± 0.15) Exp[-i(151 ± 4)°] x 10~ 6 GeV, 

E = (1.67±0.13)Exp[i(115±5)°] x 10" 6 GeV, 

A = (1.05 ± 0.52) Exp[-i(65 ± 30)°] x 10~ 6 GeV. (4.1) 

Hence, the weak annihilation (H^-exchange E or H^-annihilation A) amplitude has a sizable 
magnitude comparable to the color-suppressed internal H^-emission amplitude C with a 
large phase relative to the tree amplitude T. As discussed in [5], it receives long-distance 
contributions from nearby resonance via inelastic final-state interactions from the leading 
tree or color-suppressed amplitude. The effects of resonance-induced FSIs can be described 
in a model independent manner and are governed by the masses and decay widths of the 
nearby resonances. Weak annihilation topologies in D — > PP decays are dominated by 
nearby scalar resonances via final-state resacttering. The relative phase between the tree 
and color-suppressed amplitudes arises from the final-state rescattering via quark exchange. 
This can be evaluated by considering the t-channel chiral-loop contribution or by applying 
the Regge pole method (for details, see [5]). 

For D — > SP decays, there are several new features. First, one can have two different 
external VF-emission and internal VF-emission diagrams, depending on whether the emission 
particle is a scalar meson or a pseudoscalar one. We thus denote the prime amplitudes T' and 
C for the case when the scalar meson is an emitted particle. The quark-diagram amplitudes 
for various D — > SP decays are listed in Table IV. Second, because of the smallness of the 
decay constant of the scalar meson as discussed before, it is expected that \T'\ <C \T\ and 
\C'\ <C \C\. A noticeable example is the decay D° — > K~Oq. Its branching ratio is naively 
predicted to be of order 10~ 6 , which is strongly suppressed compared to the counterpart decay 
D° — > K~7i + in the PP sector. Experimentally, K~a^ has a branching ratio comparable 
to K~7i + . This implies the importance of the H^-exchange term in D° — > K~Oq. Third, 
since Kq and the light scalars a, k, fo, do fall into two different SU(3) flavor nonets, one 
cannot apply SU(3) symmetry to relate the topological amplitudes in D + — > f n + to, for 
example, those in D + — > K n + . Note that in flavor SU(3) limit, the primed amplitudes T' 
and C diminish in the factorization approach due to the vanishing decay constants of scalar 
mesons. 

Just as in D — > PP decays, the topological weak annihilation amplitudes E and A, 
which are naively expected to be helicity suppressed, can receive large long-distance final- 
state interaction contributions. For example, there is a contribution to the H^-exchange 
amplitude E of D° — > K* 7r° from the color-allowed decay D° — > Kq~tt + followed by a 
resonant-like rescattering. As discussed in [5] , Fig. 1 manifested at the hadron level receives a 
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X + d 




FIG. 2. Contributions to D° — ► K 7r° from the color-allowed weak decay D° — ► Kq tt + fol- 
lowed by a resonant-like rescattering. This has the same topology as the VF-exchange graph. 

s-channel resonant contribution from, for example, the 0~ resonance i^(1830) and a t-channel 
contribution with one-particle exchange. Likewise, the ^-exchange term in D° — > a^7r T 
receives the 7r(1800) resonance contribution. 



V. GENERALIZED FACTORIZATION AND ANALYSIS 



We will study D — > SP decays within the framework of generalized factorization in which 
the hadronic decay amplitude is expressed in terms of factorizable contributions multiplied 
by the universal (i.e. process independent) effective parameters a« that are renormalization 
scale and scheme independent. More precisely, the weak Hamiltonian has the form 

G 

H eS = ~^V cqi V* q2 [ai(wg 2 )(gic) + a 2 {q 1 q 2 ){uc)\ + h.c, (5.1) 

with (gig 2 ) = <?i7^(l — 75)92- For hadronic charm decays, we shall use a>i = 1.15 and 
a 2 = —0.55 . The parameters a\ and a 2 are related to the Wilson coefficients via 

at = ci(/i) + c 2 (/i) (J^- + Xi(v)j > «2 = c 2 (//) + ci(/i) (J^- + X2(/i)) , (5.2) 

where the nonfactorizable terms W1U compensate the scale and scheme dependence of 

Wilson coefficients Ci(n) to render aj physical. 



A. D K$(U30)tt 



K* -n 



K*qTi and 



Among the four measured D — > KqTI modes: D + — > K n + , D° 
— > Kq7t + , only the first one does not involve weak annihilation and hence it can be used 
to fix the form factor for D — > Kq transition. More precisely, 



A(D A 



K 7T 



^Lv V* 
' "' 



cs v ud 
2 



aiU{m 2 D - m 2 K *)F { 



DK, 



+ a 2 f K *(m D 



ml)F^(m 2 



(5.3) 
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TABLE IV. Quark-diagram amplitudes and branching ratios for various D — > SP decays with 
and without the long-distance weak annihilation terms induced from final-state interactions. Light 
scalar mesons a, k, ao(980) and /o(980) are described by the q 2 q 2 states, while Kq is assigned by 
qq. Experimental results are taken from Table I. 



Decay 


Amplitude 


E ■ 




#FSI 




^expt 




D+ - f n+ 


^«jV£(T + C" + 2A)/>/2 + K.V£ v^C" 3.5 x 10" 


-4 


see text 




(4.3 ± 1.1) x 10" 


-4 






2.2 x 10" 


-5 


2.2 x 10" 


-5 


(2.0 ±0.5) x 10- 


-4 


- a+K° 




1.7 x 10" 


-2 


1.7 x 10' 


-2 










1.7 x 10" 


-3 


1.7 x 10- 


-3 


(3.2 ±0.6)% 




— > CT7r + 


KdK*rf( T + C + 2A ) 


input 








(2.1 ±0.5) x 10- 


-3 


— > 




input 








(6.5 ± 1.9)% 






V CS F*,(T + C") 


input 








(1.8 ±0.3)% 






V cs V: d (C + 3E)/V2 


8.2 x 10- 


-4 


input for E 


(6.3 ± 1.2) x 10" 


-3 


a+tf- 


V cs V: d (T' + E) 


2.8 x 10- 


-6 


1.1 x 10- 


-3 


(2.2 ±0.5)% 






V cs V: d (C-E)/V2 


3.5 x 10- 


-3 


3.6 x 10- 


-3 


(7.9 ± 1.7)% 




a K+ 


V cd V: s (T + E) 


8.1 x 10" 


-5 


7.9 x 10" 


-5 


(2.1 ± 1.3) x 10- 


-3 


-► aj?r- 


V cd V: d (T' + E) 


1.7 x 10- 


-7 


6.5 x 10- 


-5 


(3.4 ±2.8) x 10- 


-3 


-»• a vr+ 


V cd V: d {T + E) 


1.3 x 10- 


-3 


1.3 x 10- 


-3 


(9.5 ±7.9) x 10- 


-4 


- K*-vr+ 


VcsV: d {T + E) 


1.3 x 10- 


-2 


1.1 x 10' 


-2 


(1.18 ±0.25)% 
















(7.0±?;|) x 10~ 3 




V cs V: d (C' + E)/V2 


3.9 x 10' 


-4 


3.7 x 10" 


-3 


(8.6±H) x 10- 3 






V cs V: d (2T + 2A)/V2 


input 








(1.8 ±0.3)% 






V CS V: S (2T + 3A)/V2 


1.2 x 10" 


-3 


1.2 x 10" 


-3 


(1.8 ±0.8) x 10" 


-3 




4.0 x 10- 


-4 


1.5 x 10" 


-3 


(7 ±4) x 10~ 3 






V cd V: d T + V cs V: s A 


1.3 x 10- 


-3 


1.1 x 10" 


-3 


(2.3 ± 1.3) x 10- 


-3 



From a fit to the E791 result B(D+ -> K vr+) = (1.8 ± 0.3)% (see Table I), we obtain 

Fq K *°(0) = 1.20 ±0.07. (5.4) 

As explained in Sec. II, the E791 analysis for D + — > K~it + it + has included the scalar 
contribution from the k resonance and found a better improved fit. If the PDG value of 
(3.7 ±0.4)% for the branching ratio of D + — > K rr + is employed [14], we will get a too large 
form factor of order 1.60 . 

Under the factorization approximation, the factorizable amplitudes of other D — > K^tt 
decays read 

aMml -m 2 K *)F DK °(ml) + a 2 (K* '7r + \(sd)\0}(0\(uc)\D°}_ , 



A(D° -> K*~n + ) = -^V CS V: C 
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A(D° -> <V°) = ^fV cs V: d a 2 [f KS (ml - m^)F ^(m^) + (^-vr+K^lO)^!^)^ ) 



G 



A{Dt - X *V + ) = -J= fll [ K^/^^ - ™ 2 K *)F C r° (ml) 

+ ^ s C«^ + l(^)|0)(0|(^)|L> s + )]. (5.5) 

The factorizable (or short-distance) weak annihilation contribution is conventionally argued 
to be helicity suppressed. 

We see from Table IV that the predictions for D° — > Kq~7t + and — > Kq°7v + are in 
agreement with experiment, whereas D° — > K^n and Df — > K*q K + are too small by an 
order of magnitude. The latter implies the importance of long-distance weak annihilation 
contributions induced from FSIs. If we assume that the relative phase and magnitude of 
the ^-exchange amplitude E and the ^-annihilation amplitude A relative to the external 
W-emission amplitude T are the same as in the case of D — > PP decays, namely, 



E/T « 0.63 e iU5 ° , A/T « 0.39 e^ 65 ° , (5.6) 

i^V and L>+ -> K V K+ are enhanced 
L>° -> K^n + and L>+ -> ir*°7r + are affected only slightly. 



we find that D° — > K 7r° and Z>+ — > lf are enhanced by an order of magnitude, while 



B. D + -f K7T+ 

As noticed in passing, although the decay D + — > /t7r + has very similar topological quark 
amplitudes as _D + — >• K n + (see Table IV), they cannot be related to each other via SU(3) 
symmetry as k and Kq belong to two different SU(3) flavor nonets. Using the decay constant 
f K = 65 MeV as estimated in Sec. III.B, we find that F Dk (0) = 0.85±0.10 from the measured 
D + — > ktt + rate. If the decay constant is negligible, then the form fcator will become smaller, 
F D/t (0) = 0.64 ± 0.10, due to the absence of a destructive contribution from C 



C. D+ -»• (T1T+ 

Neglecting PU-annihilation and taking m a = 500 MeV, the form factor F® a extracted 
from D + — > an + reads 

F Dct (0) = 0.42 ±0.05. (5.7) 

This is quite different from the fit value 0.8 ± 0.2 obtained in [58] using the E791 data for 
D + — > 7r + 7r + 7r _ [9] and the Breit-Wigner description of the a resonance. Note that in the 
conventional quark model, the a flavor wave function is given by (uu + dd)/\/2, while it 
is uudd in the four-quark picture. Therefore, in the SU(3) symmetry limit, the ratio of 
\A(D + — > o"7r + ) / A(D + — > K7r + )| 2 is \V c d/V cs \ 2 if a is made of four quarks, while it will be two 
times smaller if a is a qq state. 
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D. D -► / (980)tt, / (980)K 



Since the VU-annihilation contribution is smaller than the VU-exchange one in D — > PP 
decays [see Eq.(4.1)], we will neglect the ^-annihilation amplitude A as a first approximation 
and determine the D — > / form factor from experiment. We will choose P+ — > /o7r + or 
P+ — > foK + rather than D + — > /o-ff + to extract Fq °(0). The reason is as follows. In the 
SU(3) limit, it is expected that (see Table IV) 



B(D+ - f K+) 



ad 



v u . 



(5.8) 



It is easily seen that this relation is borne out by experiment. In contrast, the relation 



r(p+ -> f K+) _ 1 
t(d+ - /„#+) 4 



cd. 



(5.9) 



is different from the experimental ratio which is close to |V^/Ks| 2 - This implies that the 
decay rate of D + — > f K + inferred from FOCUS [11] is probably too large by a factor of 4. 
Indeed, since this mode is Cabibbo doubly suppressed, it is unlikely that its branching ratio 
is of the same order as the Cabibbo singly suppressed one D + — > /o 7r+ - At any rate, it is 
important to check this mode soon. From the decay — > fon + , we obtain 



r 



(0) = 0.52 ±0.04 



(5.10) 



and hence F D/o (0) = 0.26 ± 0.02 from Eq. (3.21). 

Since we have neglected the VU-annihilation contribution in the process of extracting the 
form factor P Ds (0), we will consistently ignore this contribution in all (D, P+) — > foil, f$K 
decays listed in Table IV. It is clear from this Table that one needs a sizable lU-exchange to 
account for D° — > foK°. One can utilize this mode to fix the amplitude E to be 



E/T « 0.40 e 



100° 



(5.11) 



where we have assumed a phase of 100° of the lU-exchange term relative to the tree ampli- 
tude. 

As for the decay D + — > /o7r + , the factorizable internal VU-emission amplitude is absent 
owing to a vanishing /o decay constant. Nevertheless, it does receive long-distance contri- 
bution via final-state rescattering, see Fig. 3. Since V CS V* S ~ ~VcdVud an< ^ the amplitude 
C is governed by a 2 , the amplitude V cd V* d C / V2 + V cs V* s y^C « -V cd V* d C / will give a 
constructive contribution and enhance slightly the decay rate of D + — > /o7r + . At any rate, 
the agreement between theory and experiment for this mode implies that the form factor for 
D + — > /o is indeed smaller than the one for P+ — > / . 

If / (980) is made from gg, the ratio of D + -> / ^ + to P+ -> / ^ + will be 



i 

- tan 2 9 



V cd 2 



(5.12) 
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(b) 



FIG. 3. (a) Long-distance contributions to the amplitude V C( iV* d C of D + — > /o7r + from the 
color-allowed weak decay D + — > foTr + followed by a rescattering via quark exchange, and (b) the 
long-distance contribution to the amplitude V CS V* S C of D + — > /o7r + from the color-allowed weak 
decay D + — ► i^"g + /C°, K + K*q followed by a rescattering via quark exchange. 

For # ~ 34° (see Sec. III. A), it is easily seen that the two-quark relation (5.12) is similar to 
(5.9) for the four-quark case. However, for 6 ~ 140° as favored by hadronic J/ip decays and 
the radiative decays 4> — ► f (980)7, /o(980) — > 77, the nn and ss components in the gg wave 
function of /o(980) have opposite signs. This means that the interference between the tree 
amplitude T and the W^-annihilation amplitude A in the decay, for example, — ► / (980)7r + 
is opposite in the 2-quark and 4-quark models. That is, if the interference is constructive in 
one of the quark models, it will be destructive in the other model, unless the relative phase 
between T and A is 90°. Therefore, whether or not one can distinguish between the qq and 
q 2 q 2 pictures for / (980) via nonleptonic D decays depends on the fo — u mixing angle and 
the magnitude and the phase of the H^-annihilation term. 

Finally we comment on the decays D° -> / (1370)iT, D + -> / (1370)7r+ and D+ -> 
/o(1370)7r+ which have been measured by ARGUS [17], E687 [16] and CLEO [10], by FOCUS 
[11] and by E791 [9], respectively, with the results 

f (4.7 ± 1.4) x 10~ 3 ARGUS,E687 
= I (5.9^;?) x lO" 3 , CLEO 
B(D + -> / (1370)7T + )i3(/o(1370) -> K + K~) = (6.2 ± 1.1) x 10" 4 FOCUS 
B(Dj -> / (1370)7T + )i3(/ (1370) -> Tr+Tr") = (3.3 ± 1.2) x 10" 3 E791 (5.13) 



B(D U ^ / (1370)^)B(/o(1370) - tt+tt- 



Since the branching fractions of /o(1370) — > 7r + 7r ,K + K are unknown, the individual 



21 



branching ratio of D decays into /o(1370) cannot be determined at present. Nevertheless, 
if /o(1370) is a nn state in nature, the decay Df — > /o(1370)7r + can only proceed through 
the topological H^-annihilation diagram. Hence, this will be the first direct evidence for a 
non-vanishing iy-annihilation amplitude A in D — > SP decays. The other modes D° — > 
fo(1370)K and D + — > / (1370)7r + proceed via the internal H^-emission diagram C and the 
external W^-emission diagram T, respectively. Taking into account the phase-space correction 
and noting that D + — > / (1370)7r + is Cabibbo singly suppressed while the other two are 
Cabibbo allowed, it is obvious that |T| > |C| > \A\, as it should be. 



E. D -► a (980)vr, a (980)K 

Because the primed amplitudes T" and C" are largely suppressed relative to the unprimed 
ones T and C owing to the smallness of the ao decay constant, it is interesting to notice that 
the neutral state o^K° is not color suppressed relative to the charged mode OqK~ , contrary 
to the case of D° — > tt + K~, 7t°K°. It is also anticipated that ag7r + ^> a^ir^, a relation 
that cannot be tested by the present preliminary data as they do not have enough statistical 
significance. 

Just as the D decays to a, k, / (980) and Kq(1430), one may use the channel D + — > o,qTT + 
to fix the form factor for D — > a transition. However, in the SU(3) limit, one has the relations 
(see Table IV) 

B{D+ -> a°7r+) _ 1 V cd 2 - a°7r+) 1 



-> K7T + ) 2 



V FJ 



Experimentally, i?i ~ \ and i?2 ~ 15. This indicates that the preliminary data of D + — > 
a[]7r + is too large by at least an order of magnitude. Therefore, we will instead use the value 
of 0.77 for F,f ao (0) (see Sec. III.C). For the ^-exchange amplitude we can apply Eq. (5.11). 

The results of calculations are shown in Table IV. Obviously all the predicted branching 
ratios for D — > ao^r, clqK (except for D° — > 7r + ) are too small by one to two orders of 
magnitude when compared with experiment. Note that D° — > a$ K + is Cabibbo doubly 
suppressed and it appears to be very unlikely that it has a large branching ratio of order 
10~ 3 . From Table IV we also see that D + — > a^K has the largest branching ratio among 
the two-body decays D — > a n(K). 

If we fit the D -> a form factor to D+ -> a[j7r + , we will get F Da °(0) = 3.4 and the 
large discrepancy between theory and experiment will be greatly improved. However, in 
the meantime we also predict that B(D + — > a^TC ) = 0.35 which is obviously too large. 
Moreover, it is impossible to achieve this abnormally large D — > a form factor in the quark 
model. 

It is possible that one has to apply the Breit-Wigner approximation for ao(980) to derive 
the branching ratios for D — > aovr, a^K from the three-body decays of charmed mesons. 
Furthermore, the fraction of ao(980) — > KK should be pinned down. It will be interesting 
to compare the experimental results with the predictions exhibited in Table IV. 
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VI. CONCLUSIONS 



The nonleptonic weak decays of charmed mesons into a scalar meson and a pseudoscalar 
meson are studied. The scalar mesons under consideration are a [or / (600)], k, / (980), 
a (980) and Kq(1430). The main conclusions are: 

1. Studies of the mass spectrum of scalar mesons and their strong as well as electromag- 
netic decays suggest that the light scalars below or near 1 GeV form an SU(3) flavor 
nonet and are predominately the q 2 q 2 states, while the scalar mesons above 1 GeV 
can be described as a qq nonet with a possible mixing with + qq and glueball states. 
Therefore, we designate q 2 q 2 to a, n, a (980), / (980) and qq to Kq. 

2. The topological quark-diagram scheme for D — > SP decays is more complicated than 
the case of D — > PP. One can have two different external W^-emission and internal W- 
emission diagrams, depending on whether the emission particle is a scalar meson or a 
pseudoscalar one. The quark-diagram amplitude for the case when the emitted particle 
is a scalar meson is largely suppressed relative to the one when the pseudoscalar meson 
is emitted. Moreover, the former amplitude vanishes in the limit of SU(3) symmetry. 

3. The charmed meson to k and Kq transition form factors are extracted from the 
Cabibbo-allowed decays D + — > ktt + , K ir + , respectively, while D + — > a and Df — > f 
ones are inferred from D + — > <jtt + and Df — > /o 7r+ ) respectively, based on the assump- 
tion of negligible W annihilation. We show that a large form factor for D Kq is 
expected. The relation F D+fi) = F Ds ^°/2 obtained in the 4-quark picture for / (980) 
leads to a prediction for D + — > f Q 7i + in agreement with experiment. Note that the 
value of the form factor F Do "(0) = 0.42 ± 0.05 obtained in this work is very different 
from the one 0.8 ± 0.2 quoted in the literature. It is pointed out that the ISGW model 
and its improved version the ISG W2 model predict too small D — > S form factors even 
at zero recoil. 

4. Except for the Cabibbo doubly suppressed decay D + — > foK + , the data of D — > 
an, / 7r, foK, KqTt can be accommodated in the generalized factorization approach. 
Sizable weak annihilation contributions induced from final-state interactions are crucial 
for understanding the data. For example, the importance of the W -exchange term is 
implied by the decays D° — > foK°, K*q71° and the W^-annihilation one by — > 
K K + . Without VF-exchange or VF-annihilation contributions, the decay rates of 
these modes will be too small by one order of magnitude. The branching ratio of 
D + — > foK + is predicted to be of order 10" 5 and should be tested soon. 

5. The predicted D — > ciqtt, a^K rates are too small by one to two orders of magnitude 
when compared with the preliminary measurements. It is pointed out that D + — > clqK° 
should have the largest branching ratio among the decays D — > a n, a K. 



23 



6. If /o(980) is a qq state in nature, it must contain not only ss but also uu and dd 
content. The fo — u mixing angles inferred from the hadronic decays J/i[) — > fo(f>/u, 
the radiative decay — * / 7 followed by / — > 77, and the strong coupling of f 
to -fTfT and irn are not quite compatible with each other. If 9 ~ 140°, then it will 
be possible to distinguish between the two-quark and four-quark pictures for / (980) 
in the decay, for example, — > /o 7r+ - 111 the SU(3) symmetry limit, the ratio of 
\A(D + — > (77r + )/A(D + — > K7r + )| 2 can be different by a factor of 2 in these two different 
pictures. 
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